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Abstract 

We investigate the primordial scalar perturbations in the thermal dissipative inflation where the 
radiation component (thermal bath) persists and the density fluctuations are thermally originated. 
The perturbation generated in this model is hybrid, i.e. it consists of both adiabatic and isocur- 
vature components. We calculate the fractional power ratio {S) and the correlation coefficient 
(cos A) between the adiabatic and the isocurvature perturbations at the commencing of the radi- 
ation regime. Since the adiabatic/isocurvature decomposition of hybrid perturbations generally is 
gauge-dependent at super-horizon scales when there is substantial energy exchange between the 
inflaton and the thermal bath, we carefully perform a proper decomposition of the perturbations. 
We find that the adiabatic and the isocurvature perturbations are correlated, even though the 
fluctuations of the radiation component is considered uncorrelated with that of the inflaton. We 
also show that both S and cos A depend mainly on the ratio between the dissipation coefficient F 
and the Hubble parameter H during inflation. The correlation is positive (cos A > 0) for strong 
dissipation cases where T/H > 0.2, and is negative for weak dissipation instances where T/H < 0.2. 
Moreover, S and cos A in this model are not independent of each other. The predicted relation 
between S and cos A is consistent with the WMAP observation. Other testable predictions are 
also discussed. 
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I. INTRODUCTION 



The recently released data of the Wilkinson Microwave Anisotropy Probe (WMAP) con- 
firmed the earlier COBE-DMR's observation about the deficiency in fluctuation power at 
the largest angular scales The amount of quadrupole and octopole modes of the CMB 

temperature fluctuations is anomalously low if compared to the prediction of the ACDM 
model. It implies that the initial density perturbations are significantly suppressed on scales 
equal to or larger than the Hubble radius. Models of structure formation with a cut-off 
power spectrum of perturbation on large scales provide a better fit to the CMB temperature 
fluctuations. The most likely cut-off wavelength derived from the WMAP data ^ actually 
is the same as that determined by the COBE-DMR 0, |^ . 

The super-horizon suppression is difficult to make compatible with models which produce 
pure adiabatic (isentropic) perturbations. However, it might be explained if the perturba- 
tions are hybrid. The different behavior of adiabatic and isocurvature (entropic) perturba- 
tions around the horizon scale can be used to construct power spectra with a super-horizon 
suppression. The WMAP data show not only a possible non-zero fraction of isocurvature 
fluctuations in the primordial density perturbations, also the correlation between the adia- 
batic and the isocurvature components 0. These results then turn into the constraints on 
the multi-component inflationary models, as the initial perturbations generated from these 
models are principally hybrid [7|. The double and multi- field models have been extensively 
studied in this context jsf. 

In this paper we will investigate the hybrid perturbations created by an inflation with 
thermal dissipation, the warm inflation scenario j9| . In the scheme of the thermal dissipative 
inflation the universe contains a scalar field and a thermal bath during the inflation era. The 
two components are coupled via the thermal dissipation. In addition to fitting the amplitude 
and the power law index of the power spectrum given by the COBE data [10|, the thermal 
dissipative inflation leads to a super-horizon suppression of the perturbations by a factor 



> 0.5 Recently, it has been found that the warm inflation of a spontaneous symmetry 
breaking potential with strong dissipation is capable of accommodating a running spectral 
index n of the primordial perturbations, and generally yields n > 1 on large scales and n < 1 
on small scales jl^. Our purpose here is to study the fractional power of the isocurvature 
perturbations, as well as the cross correlation between the adiabatic and the isocurvature 
fluctuations in the thermal dissipative inflationary model. 

In contrast to a single or a double field inflations, the evolution of the universe in the ther- 
mal dissipative inflation does not need a stage of non-thermal post-inflationary reheating. As 
long as the damping coefficent F satisfies the criterion given in US, r > {M/mvifH, where 
mpi, M, and H stand for the Planck energy, the energy scale, and the Hubble expansion 
of the inflaton respectively, the dissipation is effective enough to make the temperature of 
the radiation component increase continuously during the inflationary epoch. The universe 
would eventually enter the radiation-dominated phase when the temperature is high enough 
so that the radiation component prevails. Since the evolution of entropy only depends upon 
the thermal dissipative process during inflation, the entropic perturbations are not contam- 
inated by the entropy production in the reheating stage. Therefore, the primordial hybrid 
perturbations induced by the thermal dissipation can be calculated unambiguously. 

The dynamical background of the thermal dissipative inflation model has been investi- 
gated within the framework of quantum field theory. It has been shown that the dissipation 
may amount to the coupling of the inflaton to a large number of particle species Bill. In 
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this sense, the two-field model and the thermal dissipation model can be considered as two 
extremes among multi-component infiations. The former adds one more field to the single 
infiaton, while the later has a large number of additional fields. 

The adiabatic and the isocurvature perturbations in the thermal dissipative model have 
been estimated in Yet, these calculations are not immune from the problems 



induced by gauge issues which are crucial for thermal dissipative perturbations In 
particular when interactions between the infiaton and the thermal bath are substantial, 
the commonly used adiabatic/isocurvature decomposition is not gauge- independent on the 
ground of super-horizon. Therefore, we must take a full relativistic treatment to analyze 
the evolution of the hybrid perturbations generated in the thermal dissipative infiation. 
Moreover, the fiuctuations of the radiation component have not been carefully considered 
in previous works. Although the energy fiuctuations of the radiation component are always 
less than that of the infiaton field, they are not negligible in examining the relative phase 
between the adiabatic and the isocurvature perturbations. 

This paper is organized as follows. In §11 we introduce the thermal dissipative infiationary 
model in relativistic covariant form. The initial adiabatic-to-isocurvature ratio is given in 
§111. Sec. IV presents a full relativistic calculation on the super-horizon evolution of adiabatic 
and isocurvature perturbations. The numerical result of the spectrum of the adiabatic-to- 
isocurvature ratio is also given in §IV. We then summarize our findings in §V. The appendices 
provide the necessary details of the relativistic theory of linear perturbations. 



II. THE MODELS OF THERMAL DISSIPATIVE INFLATION 
A. The background field and radiation 

We consider a universe consisting of a scalar infiaton field and a radiation component 
which mimics the thermal bath. The total energy-momentum tensor is 

Tab = Ti^<t>)ab + Ti^r)ab, (l) 

where subscripts (</>) and (r) are respectively for the scalar field and the radiation, and the 
Latin indices run from 1 to 4. 

The energy-momentum tensor can be decomposed into fiuid quantities as jl^ 

Tab = PUaUb + phab + QaUb + Qb^a + T^ab, (2) 

in which QaU"" = Habu'^ = 0, Hab = TTba, and Ua can be any timelike vector field, i.e. u°'Ua = — 1, 
which is generally taken to be the average velocity vector. The total energy density of matter 
measured by an observer is represented by p, while p and iTab denote the isotropic and 
anisotropic pressures respectively. The quantity prescribes the energy fiux relative to u"". 
With Ua, the line element of the spacetime can be written as 

ds'^ = gabdx"'dx^ = — {uadx"')^ + habdx"'dx^, (3) 

where hab = gab + UaUf, is a projection tensor which maps points into the rest space of the 
observer u"'. Hence at the space-time point x", the observer u"" assigns to the event + fix" 
a spatial separation (habdx'^dx'^Y^'^ , and a time separation Uadx"' from him. 
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For a minimally coupled scalar field 0, the Lagrangian is 

L = (4) 

where V{(f)) is a self-interaction potential. The fluid quantities of the energy-momentum 
tensor for the (j) field are then 

PW = ^(</',a«")' + y = l{<l>,aU''f - V{<l>), (5) 

<i{<t>)a = -<t),cU''h\(t)fi, 7r(^)a6 = 0. 

Assuming the radiation component is a relativistic ideal fluid, we have 

p^r) = kT^, P{r) = (l/3)p(r), (6) 

Q{r)a = 0, T^{r)ab = 0, 

where T is temperature, k, = {n'^ /30)gcs and Qcs is the effective number of degrees of 
freedom at temperature T. As mentioned previously, the thermal dissipative inflation is a 
multi-component model with very high multiplication, the parameter g^s amounts to the 
actual number of components. For simplicity, we will absorb the factor k into T, and use 
= T"^ in what follows. 



B. Interactions between the scalar field and the thermal bath 

The total energy-momentum conservation of the system is governed by 

T^]b = 0. (7) 

The interactions between the field and the thermal bath are characterized by the force 
vectors defined as 

Q{4,)a = T(^<l>)a;b^ Q{r)a = 7'(r)a;6' (8) 

Obviously, we have 

Q{<j>)a + Q{r)a = 0. (9) 

The interaction term Q(i)a {i is for or r) can be further decomposed into 

Q{i)a = Q{i)Ua + J{i)a and J(i)a = 0. (10) 

The quantities and J(j)a are, respectively, the temporal and the spatial components of 
Q(i)a- They describe the energy and the momentum exchange between the scalar field and 
the thermal bath. Substituting (10) into (8), one has 

Q{i) = -u'^T^^^.^ J[i)c = KT^^)^.^. (11) 

The dissipation of the scalar field can be modeled by 

Q{4>)a = -Q{r)a = -^{(p,bu'')4),a, (12) 

where F can be a function of 0, and is always positive. With Eq. (10) we have 

Qw = -Qir) = -r(0X)'> (13) 

and 

J^c = - J(r)c = -r{(l),bu')(l),aK. (14) 
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III. THE INITIAL ADIABATIC AND ISOCURVATURE PERTURBATIONS 



A. The background solutions 

To find the background solutions, we consider all quantities being uniform and isotropic. 
Accordingly, the space-time of the universe assumes the fiat Friedman-Robertson- Walker 
metric, 

rfs2 ^ g^^dx^'dx^ = -dt^ + [dr^ + (dO^ + sin^ 9d(j)'^)] . (15) 
The Einstein equations of the expanding universe yield 

iir = -47rG(02 + 7P(.)), (16) 

and 

^ ^ [P« + + V{cf>)\ , (17) 

where 7 = 4/3 represents the adiabatic index of the thermal radiation, a denotes the cosmic 
scale factor, and H — a/ a is the Hubble parameter. 
The equation of motion for the scalar field is 

+ 3i/(/) + r(/)- a- W + l^,^ ((/.) = 0. (18) 

For the uniform background field 0, the term in (19) can be ignored, and we have 

0+(3i/ + r)0 + y'(0) =0. (19) 

where the / denotes d/d(j). 

The equation of motion for the radiation component (the thermal bath) is derived from 
the first law of thermodynamics as 

+ 3i/7P(,) = rc/)^. (20) 

Thus, the quantity F in (12) represents the dissipation "coefficient", which describes the 
production of radiation from the field. Furthermore, for the background solution, Eq. 
(14) gives 

J{d>)c = J(r)c = 0. (21) 

This is expected as the uniform and isotropic scalar field comoves with the radiation field, 
the net momentum exchange between them vanishes. 

When the potential energy V is dominant, i.e. V ^ -|- the Hubble parameter 
H{t) depends largely on V, and the universe undergoes an infiation. In the slow-roll regime 
V" < 3i/2, Eq. (19) yields 

V 

(t>^ . (22) 

The scalar field approximates the trajectory 

= 0oe''* (23) 

with 13 ~ -V"{Q)/{?,H + T) < H. Thus, with Eqs. (19) and (20), the behavior of the 
radiation component at the infiationary phase can be characterized by 

T'^p,,^%^-^L.i\ (24) 



5 



This solution actually can be obtained directly from (20) by considering the slow rolling 
condition ~ j^. Hence, Eq. (24) is still true even when T is 0-dependent. Numerical 
solutions to Eos. (19) and (20) show that (24) is indeed available when F is a power law 
function of |lO|. 

Thus, the temperature of the thermal bath increases mostly during the inflation era. In 
the standard inflation scenario which corresponds to the F = case, the radiation is blown 
off as exp(— 4ift). Therefore, the solution [Eq. (24)] is independent of the initial value of 
P(r.) when t is greater than 1/H. 

Equations (19) and (20) imply that the field dissipation would eventually heat up the 
universe, giving rise to the co-existence of the scalar field and a radiation component with 
temperature larger than the Hawking temperature, i.e. 

T> H. (25) 

Under this condition, the large scale reheating is unnecessary for inflation with dissipation. 
That is, the inflation regime will smoothly transfer to the radiation-dominated regime when 
T is high enough, and the radiation component becomes dominant [^. Numerical solutions 
to (19) and (20) have shown the smooth transition jlflj. This feature is critical for calculating 
the primordial entropic fluctuations because the initial perturbations will be unaffected by 
the large scale post-inflationary entropic process, such as the reheating. 



B. Initial perturbations of the scalar field 

The initial perturbations of the field is calculated by the linearly perturbed field equa- 
tion upon the space-time background (15). Because the primordial fluctuations are produced 
well within the Hubble radius, one can use the calculation without considering the gravita- 
tional gauge jl3|- For instance, the Fourier mode of the field perturbations with comoving 
wavenumber k is characterized by a Langevin-like equation 

d6<j> ea-^ + v"{^) 

-dr = — 3H+T ^'^+^- (2^) 

The noise term rj given by thermal fluctuations is Gaussian. The statistical property of the 
noise t] can be determined by means of the fluctuation-dissipation relation |9| 

(V) = 0, (27) 
{r]{k,t)r]{k',t')) = D5k,-k'5(t - t'), (28) 

with D = {3H^T/2'k){3H + r)~^, and (...) denotes averaging over an ensemble. The relation 
that D = {H^/2tt) for the quantum case without dissipation is easily recovered. 

Taking the slow-roll condition \V (0)| <^ 9H^ into account, the V (0) term in Eq. (26) 
can be ignored at the horizon-crossing where a/k = H~^. Accordingly, the correlation 
function of the fluctuations is given by 

(50(t)50(t'))-/^^^e-(-*')^^/(3^+n, t>t\ (29) 
and the amplitude for the horizon sized perturbation 50 is about 

-Ht) . (30) 
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which imphes the thermal fluctuations dominate over the quantum ones in the era of T > if. 
Taking the ensemble average, we obtain the perturbation 50 just outside the horizon as 



-M- 



3H + r 

The perturbations in the fleld energy density is then 

(5p(<^) = 050 + V'6(f) - V'6(f) = -{3H + r)0(50 ~ - [2(ri^ + 3)TJj^^^ 
where we have used the slow-roll condition (22), and Tfj = T/H. 



(31) 



(32) 



C. Initial perturbations of the thermal bath 

We consider the radiation component in thermal equilibrium with a temperature T . At 
sub- horizon scales, fluctuations in the thermal radiation can be estimated by 5p(r)/ P(r) — 
1 / ^Jn^ , where denotes the total number of the relativistic particles within the horizon 
H^^. Since the photon number density is proportional to T^, and the volume within the 
Hubble radius is about (47r/3)i/~^, we flnd that 



Pir) 

Accordingly, the energy fluctuations caused by 5p(>r) are characterized by 




(33) 



Sp ^ ^ J ±HT ■ HT^ = HT^Scj), (34) 

P(r) V 47r 

where Eq.(30) has been used. It should be emphasized that (34) qualifles only the relation 
between the variances of 5p(r) and 6(f), but not their phases. Principally thermal fluctua- 
tions of radiation component is random-phased with respect to the fleld fluctuations, i.e. 

(505p(r)) = 0. 

With the help of the background solutions (22) and (24), the relation between the fluc- 
tuations in the thermal energy and the fleld energy can be established, 

HT^ pi/2 

%0 = Vw^ = ^(^;^^PW- (35) 

Therefore, Sp(^r) is always less than 5p(^) in situations either T < H ot T > H. Apparently, 
the perturbed cosmic matter is just about isothermal. 



D. The adiabatic vs. the isocurvature initial conditions 

The energy perturbations of the fleld and the radiation component can be decomposed 
into adiabatic (ad) and isocurvature (en) modes as 

Spw = ^Pw+M^), (36) 
5 Pir) = Spp^+dp'^^y (37) 
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By definition, the adiabatic perturbations 5^ are given by 



P{4>) + P{<t>) P(r) + P{r) ' 

while the isocurvature mode satisfies 

Therefore, the adiabatic perturbations can be rewritten as 



rad 



Spii + Spir) 



+ h 



(r) 



PW + + P{r) + P{r) 

and the entropic perturbations are characterized by 

^PW ^Pjr) _ Sp{^) Sp^. 



PW + PW P{r) + P{r) h 



{</>) 



h 



(r) 



(38) 



(39) 



(40) 



(41) 



where we have used /i(<^) = = 0^ from Eq.(5), and h{r) = p[r) + P{r) — ^h^^/^ 

from (6) and (24). The radiation part Sp(^r)/h(r) in the above definition (41) is noteworthy 
as it may be compatible to the inflaton part 6p(^)/h(^), even though the condition 5p(r) < 
is always fulfilled. That is, 6p(^r) must not be ignored when treating the isocurvature 
perturbation no matter how small it may be. 

Considering the case when (5p(^)5p(r)) = 0, one can obtain the variances for the adiabatic 
and the entropic perturbations by virtue of Eqs. (30), (32) and (35): 



[(((5^d)2^]l/2 ^ 



1+ 



H 



.1/2 



4(r^^ + 3)2j 



PW 



h 



1 + 



9r 



1/2 



H 



4(r^, + 3)^ 



The correlation between the two perturbation components is given by 

/ \ 2 

^^en^ad\ 



. h't') I 



4(r^ + 3)2 



1 + ^ 

3 



(42) 
(43) 

(44) 



Equation (44) shows that even when the fluctuations of the held and the radiation are 

uncorrelated, the correlation between the adiabatic and the entropic perturbations can be 
significant. In these initial conditions all and 0-dependent quantities as well as V , h(^^^, 
and Fff are taken to be their values at the /c-mode horizon- crossing times. 



IV. SUPER-HORIZON EVOLUTION OF ADIABATIC AND ISOCURVATURE 
MODES 

A. Adiabatic/isocurvature decompositions on super-horizon scale 

In super-horizon regions, physical quantities should be expressed in a gauge-independent 
fashion. However, the adiabatic/isocurvature decomposition of perturbations shown in the 
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last section are not gauge-independent, i.e. either eqs.(40) or (41) are not gauge invariant 



variables. As usual |20|, |2l|, we choose the condition to fix the coordinates and to make (40) 
and (41) physically meaningful. As shown in Appendix C, a gauge invariant (GI) variable 
C can be defined as 



c = <^ + 



(45) 



3(ptot + Ptot) Hh<p) + hr)) ' 

where ip is the perturbed variable of 3-space curvature [see Eq. (A2)]. Thus, if we fix 
the coordinates by = 0, the so-called uniform-curvature gauge (UCG) condition, then ( 
becomes 6^*^ [Eq.(40)]. Thus, the initial condition and the evolution of the super-horizon- 
sized adiabatic perturbations can be described properly by the equations of ( in the UCG. 
On the other hand, we can define another GI variable w using Eqs.(Cll) and (C12) as 



6J, 



h 



(0) 



h 



(r) 



(46) 



where SJ(^) and 5J(r) is the linear perturbed variables of energy exchange between the (f) 
field and the thermal bath [Eqs.(14) and (A17)]. Evidently, w is nothing but 5*^° if replacing 
(5p(^) and Sp(^r) in (41) respectively by 5p(<^) + and 6p(^r) + ^^(r)- We may call zu a 
modified isocurvature (entropic) perturbation. As a consequence, there exists no suitable 
coordinate fixing to make both 6'^'^ and S'^'^ simultaneously a GI variable as there is an 
energy exchange between the background components. In this regard, the decomposition of 
perturbations by (40) and (41) at super-horizon scales does not have clear physical meanings. 
Instead, one should decompose the super Hubble perturbations into the adiabatic and the 
modified isocurvature modes. The initial condition and the evolution of the modified entropic 
perturbation outside the horizon can then be described unambiguously by the equations of 
uj in the UCG. Moreover, at the commencing of radiation regime when the infiation ends, 
0^0, the energy ex change in the background ceases and (46) reduces to (41). Therefore, 
perturbations defined by Eq. (46) can be explained as the isocurvature perturbations at the 
onset of the radiation-dominated epoch. 

The initial condition of 5J(0) can be estimated via (A17) as 



(47) 



Comparing to (30) and (32), 5J(</,) is always less than 5p((/,), but it can be larger than 6p(^r) if 
T > H. Since + SJ(^r) = 0, we have 5p(</,) + 5</(0) + ^P(r) + SJ{r) = ^Pi^j}) + ^P{r)- Hence, 
(40) remain vahd even when 5p(0) and Sp^r) are replaced by 5p(0) + and Sp(^r) + 5J{r)- 
More importantly, since 5 J(r) involves the effect of 0, it is correlated with the (f) field fiuc- 
tuations. The perturbations 5J(j.) and 5p(0) actually are in phase. Therefore, it is interesting 
to consider the case of imposing the phase correlation into the initial conditions. Connecting 
and 5p{r) by the relation (35), Eqs. (40) and (41) yield 



c-ad 



h 



1 + 



.1/2 
H 



2(rH + 3) 



1 + ^ 

3 



h 



(</-) 



3r 



-1/2 



H 



2{Th + 3) 



(48) 



(49) 
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and 

Consequently, 5*^" and 6'^'^ are in phase when Th > 0.217, but are anti-correlated when 
< 0.217. 



3r 



'1/2 



H 



2(rH + 3) 



ri/2 

^ H 

2(r^^ + 3) 



B. Equations of the super-horizon evolution of perturbations 

The evolution of the linear perturbations outside the Hubble radius is described by a set 
of equations of all perturbed matter variables, such as 5p(0) and (5p(r), as well as perturbed 
variables of the space-time metric. The sophisticated formalism of deriving these equations 
are given in Appendices. For the UCG, these equations are shown as (Dl) - (D8) in Appendix 
D. With the solutions of Eqs. (Dl) - (D8), it is straightforward to trace the evolution of 
relevant variables of the perturbations. 

The evolutionary features of the super-horizon-sized perturbations can be seen from Eqs. 
(D7) and (D8), which are 



50 + 



(3if + r) - 

327rG^ 



SttG, 



H 



3 H 



50 + 

-2V' 



3 H 



[X + Hx] 



(51) 



X + 



H + 



-inG 



327iGV 



3 H 



X + 



A-rG ( -V 



2^ 
3^ 



3H ^ 3H 



X 



(52) 



These equations portray the linear evolution of the A;-mode perturbation of the scalar field 
(50, and that of the shear x of the space-time metric [20]. 

The Friedmann Eq. (17) implies (3if + r) ^ {^tiG /3){(fP' / H). For modes at sub-horizon 
scales (A;Va2) > H"^, we have \{k'^la^) + > (167rG/3)(01/7i^)]. Accordingly, if taking 
the right hand side to be zero, Eq. (51) within the horizon is exactly the same as the 
linearized Eq. (18), or (26) with the slow-rolling condition. Since the initial perturbation 
of field is sub-horizon-scaled, and is governed by (18) or (26), it is consistently to assume 
the initial conditions for x and x to be zero. 

Equations (51) and (52) can be regarded as two coupled oscillators with a time dependent 
mass, damping coefficients and coupling coefficients. Whenever the "mass" becoming nega- 
tive, the perturbations undergo a decaying or growing process. This gives rise to the grav- 
itational clustering. The "mass" of the x-oscillator generally is positive for modes beyond 
the Hubble radius {k"^ /o?) <^ H^. Therefore, the perturbations 50 will not be magnified by 
gravity outside the horizon, and will retain approximately their initial values. Consequently, 
variations in both 5^ and 5'^^ are insignificant during their super-horizon journey. 



10 



C. Numerics 



We calculate the power fraction S of isocurvature perturbations, and the cross-correlation 
cos A between the adiabatic and the isocurvature perturbations at the end of the inflation 
epoch when the universe enters the radiation regime. These two quantities are defined by, 
respectively, 

1 ,,,, 

cosA = .^^^^ ' (54) 

Two sets of initial conditions, (42) - (44) and (48) - (50), are imposed respectively to 
calculate S and cos A. Both quantities depend mainly upon the dissipation parameter 
Fff. The results are plotted in Figs. 1 and 2. Figure 1 shows the prominent mixture and 
correlation of two distinguished perturbation components for the case where the fluctuations 
of the field and of the thermal bath are totally uncorrelated. The F/^-dependence of the 
curves for S and cos A in Fig. 2 reveals the similar behavior as in Fig. 1, but with different 
amounts. In particular, the isocurvature perturbations and the adiabatic ones are in phase 
when Th > 0.217, but are anti-correlated when Th < 0.217. 




0.0 



0.01 



0.10 



1.00 



10.00 



FIG. 1: The power fraction S of the isocurvature perturbations, and the cross-correlation cos A 
between 6^*^ and 6'^'^ as a function of parameter Th at the end of thermal dissipative inflation for 
initial conditions (42) - (44). 
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FIG. 2: Same as Fig. 1 with initial conditions (48) - (50), i.e. taking into account the correlation 
given by the energy exchange between the background components. 

V. CONCLUSIONS AND DISCUSSIONS 

We show that the thermal dissipative inflation produces the initial hybrid perturbations 
with signiflcant correlation. This is largely due to the coexistence of two components, 
the radiation and the (j) fleld, during the inflationary epoch. The evolution of the density 
perturbation of radiation and that of the (p held on scales larger than horizon are governed 
by different equations (D3) and (D7). Consequently after reentering the horizon, the density 
perturbation of radiation and that originated from fleld generally are different. The super- 
horizon analysis is essential to reveral the formation of the hybrid perturbations. 

We have calculated the power fraction of the isocurvature perturbations (S), and the 
correlation between the adiabatic and the isocurvature perturbations (cos A) at the end 
of the inflationary epoch. Since the transition from the inflationary era to the radiation- 
dominated period is smooth without an intervening reheating process, the values of S and 
cos A can be directly used as the initial conditions for the radiation regime. 

We found that, for each k-mode, S and cos A are mainly determined by the parameter Th 
at the instant as the perturbation mode crosses outside the Hubble horizon. That is, S and 
cos A depends entirely upon the dissipation of the inflaton. When taking into account the 
effect of energy exchange, cos A is more sensitive to the dissipation parameter Th- For strong 
dissipation cases where F > 0.217// the adiabatic and the isocurvature perturbations are in 
phase. Under weak dissipation F < 0.217if, however, the two perturbation components are 
anti-correlated. 
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Given that the current observational constraints on 5* and cos A are rather diverse |6l. l22L 
2^, we do not make the detail parameter fitting in this paper, but disscuss some properties 



of the hybrid perturbations which are useful for further model testing. Apparently, in this 
thermal dissipative model S and cos A are not independent of each other. The predicted 
S-cosA relations can directly be seen from Figs. 1 and 2. These relations hold for all k- 
modes, but they are not associated with the parameter F. On large scales these predictions 
can be confronted with the CMB observation without considering the effects of evolution 
of perturbation during the radiation era. The relation cos A ~ ±\/S shown in Fig. 2 can 
then be tested by the observed adiabatic/isocurvature ratio and correlation. For instance, 
to realize the thermal dissipative inflation, the dissipation parameter Th can be taken in the 
range ~ 0.13 — 0.4 [1, Il3,[l3- Consequently, S amounts to about 10% and cos A ~ —0.3 to 
+0.3. 

Secondly, the ^-dependence of cos A in this model is governed by the fc-dependence of 
Th- On the other hand, Eqs. (42)-(43) or (48)-(49) shows that the two differences 

dln{{6^y) rfln((,5"")^) £ \n{{6'^y) £ In ((,5"")^) 

dink dink dlnk"^ dlnk'^ 

are also specified by the fc-dependence of Th- Therefore, the difference between the spectral 
indices, or the running spectral index, of the adiabatic and the isocurvature perturbations 
should also be determined by the /c-dependence of the correlation cos A. For instance, if 
cos A are /c-independent for some range of k, the spectral indices should be fixed without any 
changes in that range. On the contrary, if there exists difference between the spectral indices 
of the two perturbation components, we should see a fc-dependent cos A. This property is 
a robust prediction regardless of the initial conditions [Eqs. (42)-(44) or (48)-(50)] in use, 
and is effective to testing the thermal dissipative inflation model. 
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APPENDIX A: PERTURBATION VARIABLES 

Appendixes A - D will give all formulas needed for a full relativistic theory of linear 
perturbations for a system consisting of a field and a fluid, such as radiation. Most 
material is taken from [2fl, . 



1. Perturbing the space-time metric 

The perturbed metric with respect to (jab [Eq.(3)] is gab = Qab + ^Qab- In this paper, we 
are interested only in the scal ar-ty pe perturbations. For this purpose, the perturbations can 
be described by four variables [20| 

lapse function a : goo = — a^(l + 2a) (Al) 
3 — space curvature (f : TZ = (A2) 
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expansion scalar k : 6 — 3H — n (A3) 

|7 
|7 



shear x: = X|a/3 - l^ajx'u (A4) 



where the comoving spatial metric tensor g^^^ is defined by Qajj — o^g'^^p. The vertical bar 
indicates a covariant derivative based on g)^^l . 

2. Perturbing the scalar field 

Prom Eq.(5), the linearly perturbed energy density is 

5p^^)^U{ct>,au''f + V'5ct>. (A5) 
The first term on the r.h.s. of (A5) can be calculated by 

u^) = + 0,„ 5u'' = 50,0 u' + <t>,^u'{ — \. (A6) 

By means of (Al), we have 5vP /vP ~ Sgoo/'^goo — —oi- Therefore, (A6) yields 

S{(l),au'') = S<i) - ^a. (A7) 

Thus, (A5) becomes 

5p(0) = (j)d(j) - (p^a + V,^ 5(1). (A8) 
Similarly, the perturbation in the pressure of the field is 

^P{<j>) = 0(^0 - - V,^ (A9) 

For the background solution, q(^)a = 0. Since q(^)aU°' = 0, the linearly perturbed q(ip)a = 
—(t>,c U''^h'^<P,b can be described by 

Sq{i,)a = V'W.a, (AlO) 

where 

^^(0) = -050. (All) 
Since 7r(0)a(, = [Eq. (5)], the perturbed anisotropic pressure S-Ki^^-^^ah — 0- 

3. Perturbing the thermal bath 

From Eq. (6), the perturbed variables of the thermal bath are given by 

perturbation of the energy density : e{r) = 5p(r) (A12) 

perturbation of isotropic pressure : 5p{r) — \5p(r) — \^{r) (A13) 

energy fiux : 5q(r)a = V'(r-),a (A14) 

anisotropic pressure : 5T:{r)ah — 0. (A15) 

where '?/^(r),a is the energy density fiux of radiation. 
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4. Perturbing the interaction terms 

Using Eqs. (13) and (A6) the perturbed temporal component of the interaction term is 
5Q(4>) = -^Q{r) = -V5{(t>,au''f = 2r(0(50 - ^''a). (A16) 
Similarly, the perturbed spatial component can be described by 

5J(^) = -5J(^) = r050. (A17) 

APPENDIX B: EQUATIONS OF PERTURBATION VARIABLES 

Based on the perturbed variables given in §A, the evolution of these variables is governed 
by the following equations: 
Definition of k 

K = -3v? + 3i/a; + — X, (Bl) 

(X 

ADM energy constraint 

HK-—ip^ -inG (s^r) + (pScP - (f'a + V,^ , (B2) 

momentum constraint ^ 

K-^X^ ISttG U^r) - M , (B3) 
a 

ADM propagation 

X + Hx^a + if, (B4) 

Raychaudhuri equation 

k + 2Hk ^{^- oi + At^G \2e^r) + 4050 - 40^0; - 2V,^ 50] , (B5) 

energy conservation of radiation 

e^r) + 4.He^r) = ^V'(r) + (4/3) p(^) {k - 3Ha) + 2r050 - Tcfi^a, (B6) 
momentum conservation of radiation 

iP^r) + ^HiP{r) = -(4/3)p(,)a - (l/3)£(,) - r050, (B7) 
energy conservation of scalar field 

50 + (3H + r)50 +(^^ + V"^ H = + a) - [(3// + r)0 + 2V'] a, (B8) 
The momentum conservation of scalar field yields an identity. 
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APPENDIX C: GAUGE INVARIANT VARIABLES 



Considering a gauge transformation x° = x'* + ^", in which is an infinitesimal quantity, 
any tensor quantity T transforms as T = T — C^T, where is the Lie derivative of T in 
the 4-vector field We may split the tensor T into its background and perturbed values; 
i.e., T = T + (5T. The transformation of the perturbed parameter becomes 6T = 6T — C^T. 
Applying the transformation to the metric tensor, we have 

Sgab - 5gab = -^^9ab = -(^a;6 + 6;a) = -i9acC,b + 9cbC,a + 9ab,cC) ' (^1) 

Since the background 3-space is homogeneous and isotropic, all the perturbation vari- 
ables are gauge independent under purely spatial gauge transformations. The perturbation 
equations presented above are also independent of the spatial gauge transformation. Under 
the temporal transformation factor = a^°, the perturbed variables are given by 2^, 2ll |: 



a = a-f^ , (C2) 
if = cp-HT^, (C3) 
X = X-T^, (C4) 

K = + I^Sif - , (C5) 

Spi = hi - PiT^, , (C6) 

i^i = + {p^ + Pi)T^ (C7) 

SQi = 6Qi - QiT^ (C8) 

SJ, = 6 J, + Q,T^. (C9) 

Equations (C6) to (C9) are applicable to components (0), (r) as well as the total fiuids. 

Based on the time-gauge transformations [Eqs. (C2)-(C9)] and Background equations of 
motion [Eqs. (19) and (20)], one can construct the following gauge invariant variables: 

C-^+ ,, , =^+ ,, ln*>+JPir, (CIO) 

3(Ptot + Ptot) 3(p(<^) + + p^r) + V(r)) 

3(p(r-) +P(r)) 

= ^ ^ (C13) 

APPENDIX D: EQUATIONS OF LINEAR PERTURBATION IN THE UCG 

Using the uniform- curvature gauge (UCG), the equations of the perturbed variables can 
be obtained from Eqs. (Bl) to (B8) by setting up = 0. We have 

a = x + Hx, (Dl) 
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K = 3Hx + \3H^ + 



X 



3IP 
4^ 



jx + 



3if2 + _ 



H 



4 
3 



X + 



4 A;^ ... 



4 1 



{3H + r) - 



3 H 



327iG 



vj- + (r0 + 2V' 



X , 



X 



r„ 32irGVl 




47rG' ( 




X- 







2e 
3^ 



X 



-47rG 



3/7 ^ 3/7 



From (D1)-(D4), we have 

47rGa2 r 



X 



and 



H 



X 



47rG 



(D3) 

(D4) 
(D5) 
(D6) 



(D7) 



(D8) 



(D9) 
(DIO) 
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